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An effective energy tensor for gravitational radiation is identified for uniformly expanding flows 
of the Hawking mass-energy. It appears in an energy conservation law expressing the change in 
mass due to the energy densities of matter and gravitational radiation, with respect to a Killing-like 
vector encoding a preferred flow of time outside a black hole. In a spin-coefficient formulation, the 
components of the effective energy tensor can be understood as the energy densities of ingoing and 
outgoing, transverse and longitudinal gravitational radiation. By anchoring the flow to the trapping 
horizon of a black hole in a given sequence of spatial hypersurfaces, there is a locally unique flow 
and a measure of gravitational radiation in the strong-field regime. 
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Introduction. Gravitational radiation is a prediction of 
Einstein gravity, whereby relatively accelerating masses 
produce ripples in space-time which propagate outwards 
with the speed of light. The waves predicted to reach 
Earth from distant astrophysical sources are so weak that 
they have not yet been directly observed, though strong 
indirect evidence comes from observations of the Hulsc- 
Taylor pulsar, for which the orbital period is decaying at 
the predicted rate due to loss of energy by gravitational 
radiation [1|. To make direct observations, a new gener- 
ation of gravitational-wave detectors is being developed, 
promising to open a new window on the universe. 

A theoretical problem has been the difficulty to define 
gravitational radiation in general. In the short-wave ap- 
proximation for linearized gravitational waves, one sep- 
arates space-time ripples from background curvature by 
assuming that appropriate length scales exist, leading to 
an averaged effective energy tensor 0. At infinity in an 
asymptotically flat space-time, an exact characterization 
was found by Bondi and others |3j, with the energy flux 
of gravitational radiation contributing to an energy-loss 
equation, whereby the total mass or energy of the space- 
time decreases due to outgoing radiation. Despite these 
successes, they reflect weak-field definitions. In partic- 
ular, since one of the promising sources is the inspiral 
and merger of a binary system of black holes or neu- 
tron stars, there has been an intense effort to predict 
waveforms using numerical simulations. Such simulations 
cover the strong-field regime of General Relativity, while 
waveforms have had to be extracted subsequently. 

This Letter reports an effective energy tensor for grav- 
itational radiation which can be applied in the strong- 
field regime, entering an energy conservation law in the 
same way as the matter energy tensor. At infinity, they 
respectively recover the Bondi energy flux and energy- 
loss equation, while on a black-hole horizon, they have 
the same form as the effective energy tensor and energy 
conservation equation recently found for dynamical black 
holes |4L the latter originally derived by Ashtekar & Kr- 
ishnan jfjj . The measure of gravitational mass-energy en- 
closed by a spatial surface is taken to be the Hawking 



energy |(|, which was useful in both cases above, reduc- 
ing respectively to the Bondi mass and irreducible mass. 
Attempts to find similar results for other definitions of 
energy have not led to justifiable physical interpretation. 

The main intended scenario is one of the above simu- 
lations, where a family of spatial hypersurfaces has been 
numerically evolved in time, and one or more marginal 
surfaces have been found in each hypersurface, practi- 
cally locating the black hole(s). A marginal surface is to 
be propagated outwards in the hypersurface, generating a 
one-parameter family of surfaces, such that the expansion 
(of the surfaces in the hypersurface) is constant on each 
surface. It transpires that the dual flow, explained below, 
should also have constant expansion. These flows were 
previously described as uniformly expanding Q and gen- 
eralize inverse mean-curvature flows |8j . Then the energy 
is found to be monotonically increasing along the flow, as 
described recently by Mars 0, and to satisfy an energy 
conservation equation as described. Further details and 
developments are to be presented elsewhere [T(i| . 

Spin- coefficient formalism. It is convenient to use the 
spin-coefficient formalism [Tlj . since it is easier to un- 
derstand the interpretation of the components as energy 
densities of ingoing and outgoing, transverse and longi- 
tudinal gravitational radiation, by comparing with the 
Szekeres characterization 0] m terms of the Weyl spinor. 
Since the spin-coefficient formalism is not generally famil- 
iar, a brief summary is given below, referring to Penrose 
& Rindler jll| for details. Henceforth equation numbers 
in triples will indicate equations in that reference. It is 
necessary to generalize some expressions and convenient 
to modify some notation. 

The basic geometrical object is a spin-metric, an an- 
tisymmetric 2-form sab (2.5.2) acting on complex 2- 
spinors (spin- vectors). It can be regarded as a square 
root of the space-time metric, which is expressed as 
fjab = £ab£A'B', (3.1.9), where the bar denotes the com- 
plex conjugate and the abstract indices indicate how a 
tensor can be expressed as a spinor dyad. A spin-basis 
is an ordered pair of non-parallel spin- vectors (o A ,i A ), 
which defines a null tetrad, a vector basis (3.1.14): 
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(l a ,m a ,m a ,l' a ) = {o o ,(ri ,60 ,i 1 ). Then 
(/,/') are real null vectors, whereas 771 is a complex vec- 
tor encoding transverse spatial vectors. In terms of the 
complex normalization factor x — £abo a l b (2.5.46), one 
finds g(l, V) = —g(m, fh) = XX with the other eight inde- 
pendent (symmetric) contractions of the tetrad vectors 
vanishing. Here the metric convention is (H ). 

The inverse metric is g^ 1 = 2(1 ®V — m®fh)/xx, where 
<g> denotes the symmetric tensor product. The dual ba- 
sis of 1-forms (4.13.32) will be denoted by (n, w, w, n') = 
g(l',—fh,—m,l)/xX) so that n(l) = w(m) = w(fh) = 
n'(l') = 1, with the other twelve such contractions 
vanishing. Then the metric is g — 2xx( n ® n ' ~ 
w <8> w). The tetrad covariant derivative operators are 
denoted by (D,5,S',D') = (V/, V m , V fl , Vj<) (4.5.23). 
The complex spin-coefficients (k, a, p, r, e, /?, a, 7) and 
(k! , a' , p' ,t' ,s' , 0' , a' ,7') are defined by (4.5.21) and en- 
code the Ricci rotation coefficients. Further defini- 
tions are of weighted derivative operators (P, 5, 3', P') 
(4.12.15), a spinor QabA'B' and a scalar A encoding 
the Ricci or Einstein tensor (4.6.24), the Weyl spinor 
^ aba' B' (4.6.35) encoding the Weyl tensor, and its five 
complex components ^0, . . . $4 (4.11.9). 

Uniformly expanding flows. It is convenient to choose 
null coordinates (x,x') such that (n,n') — (dx,dx'), 
which implies the dual-null gauge conditions [pj = 
k = p — p = e + e = T — a — = k' = p' — p' — e' + e' — 
t' — a! — /3'. Then the formalism describes two families 
of null hypersurfaces labelled by x and x' , intersecting 
in a two-parameter family of transverse spatial surfaces. 
The main application is to a spatial hypersurface foliated 
by transverse surfaces, which generates a locally unique 
dual-null foliation, intersecting in the transverse surfaces. 

Consider a one-parameter family of transverse surfaces 
with generating vector £, normal to the transverse sur- 
faces. Then £ = £ I + £}V , where the components £ a 
are constant on the transverse surfaces. It is useful to 
introduce its dual vector £* = ££ in the normal space , 
where S is the vectorial Hodge operator of the normal 
space, with orientation chosen so that £* = £ / — £ 
Since the expansions along (I, V) are —2(p, p'), the condi- 
tion for the flow to be uniformly expanding is that (p, p') 
are constant on the transverse surfaces. 

The flows can also be expressed in terms of the unit- 
expansion vector (inverse mean-curvature vector) [l4| 

r) = -H/g(H,H) H = -g-\d\og(*l)) (1) 

and its dual 77* = *n, where H is the expansion vector 
(mean- curvature vector) [l5T | , the signs reflect the metric 
convention and *1 denotes the area form of a transverse 
surface. Then H = 2(p'l + pl')/xX, V = -(l/p+ ''/V)/ 4 , 
7/* = —(l/p— V /p')/A and the expansions along (77,77*) 
are found to be (1,0). Thus £ = a(n + cq*) generate 
all uniformly expanding flows, where (a, c) are constant 
on the transverse surfaces . This freedom suggests that 
such a flow should locally exist for any untrapped or mean 
convex surface 0, pp' < 0, in any hypersurface, though 
no proof exists. 



Energy propagation. It is convenient to use the area 
radius 



R = \/ Aj Am A = <p *1 



(2) 



where A is the area of a transverse surface. The Hawking 
quasi-local energy can be written as 



E =£;( 1 + ikt ' ! » !J " 



where G is Newton's constant. Since g(H, H) 
it can be generalized to 



E = 



R 



K + pp 1 



XX 



where 



K = 00' - pp - * 2 + $0 



n 



(3) 

&pp'/xx, 

(4) 
(5) 



is Penrose's complex curvature (4.14.20), such that the 
Gaussian curvature is (K + K)/xx, and n = XX^- 
(4.11.7). Penrose's complex generalization (4.14.42-43) 
of the Gauss-Bonnet theorem is 



K 

XX 



K 

* — - = 2n(l — 7) 

XX 



(6) 



where 7 is the genus of the transverse surface. 

The propagation equations for the Hawking energy can 
be found from the compacted spin-coefficient equations 
(4.12.32), using (P, P')*l = -2(*p, *p'), as 



DE = 



D'E = 



R 
A^G 



*p- 



K + pp' - t't' 



lOl'O' 



3n 



*p 



OOO'O' 



R 
4^G 



XX 

I K + pp' — TT 



XX 
*p 



K + pp' 



(7) 



*p 



XX 

$oio'i' - 3n 



11' 1' 



XX 



XX 



K + pp' 



XX 



(8) 



For a uniformly expanding flow, p and p' may be taken 
outside the surface integrals. One then sees that the 
terms in K + pp' cancel, and the terms in 5V and 9't 
integrate to zero by Penrose's version of the Gauss diver- 
gence theorem (4.14.69), so that the energy propagation 
equations simplify to 



DE 



D'E = 



R 
A^G 



R 
A^G 



p'jao- + <frooo'o') 

XX 

p(r'f' + $ioi<o' + 311) 



XX 



p(o-'<j' + <A> 1 iyy) 



XX 
p'(tt + $ 



010'!' 



3n) 



XX 



(9) 



(10) 
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Effective gravitational-radiation energy tensor. Then 
one sees that the non-matter terms can be encoded in a 
spinor H corresponding to <f> by 



— OOO'O' 

Sioi'o' : 



~J\2 



—010' 1' 
—lll'l' 



(11) 



which can be written in a more manifestly invariant form 
as 

^ABA'B' = iBE'DClB'ED'C x (12) 

{o A o A ,o C o C 'o D o D 'l E l E ' + l A Z A ,l C T?\ D Z D 'o E o E ') 

where r ) AA 'BC are the spin-coefficients (4.5.2). Then 5 
can be seen to be a spinor in the normal space, since it 
involves (o, t) only in the combinations I — oo, V = il, 
not m — oT, fh — lo. By comparing with the Einstein 
tensor G a b = —2$ABA'B' — 6As A b£a'B' (4.6.24) and the 
Einstein equation G a b + 6Ag a b = —8TrGT a b (4.6.30), the 
effective gravitational-radiation energy tensor O is given 
by 



AirGOab = zlaba> 



(13) 



The positivity of the components of © indicate that the 
gravitational radiation carries positive energy; more ex- 
actly, O satisfies the dominant energy condition. 

Introducing a non-standard but logical index corre- 
spondence a = (0,3,2,1) <-► AA' = (00', Of, 10', 11'), 
the components 6 o = \o\ 2 /A-kG and 8n = \a'\ 2 /AnG 
have the same form as the energy densities of transverse 
gravitational radiation at past and future null infinity 
respectively 0, taking the usual convention that I 
is future outward and V is future inward. On physical 
grounds, one might expect the other components to be 
longitudinal gravitational terms. As a first step, consider 
the initial data for the dual-null initial-value problem, 
based on a transverse surface S and the null hypersur- 
faces (£, £') generated from it in the (I, I') directions. 
From the form of the spin-coefficient equations (4.12.32) 
and Bianchi identities (4.12.36-41), the free vacuum ini- 
tial data can be prescribed as (p, p',^) on S, (a, t') on 
£ and (ct',t) on Then it is natural to interpret the 
last two pairs as ingoing and outgoing radiation profiles 
respectively. 

Szckeres ^3 introduced a gravitational compass, a 
tetrahedral arrangement of four test particles joined by 
springs, to characterize the physical meanings of compo- 
nents of the Weyl tensor by their effect via the geodesic 
deviation equation. In terms of the Weyl spinor ^, 
this simplifies to: for ingoing transverse radiation, 
for ingoing longitudinal radiation, ^2 for the quasi- 
Newtonian (Coulomb- like) gravitational field, ^3 for out- 
going longitudinal radiation and "J 4 for outgoing trans- 
verse radiation. A comparison can be made using the 
recent unified framework for null and spatial infinity in 
terms of advanced and retarded conformal factors (u>, u/) 
(l(| . For future null infinity, one can use a radial null co- 
ordinate r = uj^ 1 such that Pr = 1/V2, P'r = 0. Quot- 
ing the asymptotic behaviour p ~ — l/v2r, p' ~ l/y/2r, 



a ~ aor , r ~ r r , r ~ r r , a ~ a' r 1 and as- 
suming the standard "peeling" behaviour \Ev ~ ^r e ~ 5 
(9.7.38), (4.12.32c) yields r = -f^. Then the primed 
versions of (4.12.32b,c) yield P'a' Q = f a 4 and P't = 
Thus a wavelike profile of ^3 or ^4 in retarded time 
x' will drive a wavelike profile of r or a' respectively. 
With similar results at past null infinity, this corrobo- 
rates the proposal that the radiation may be encoded 
in the relevant spin-coefficient, specifically er for ingoing 
transverse radiation, r for outgoing longitudinal radia- 
tion, t' for ingoing longitudinal radiation and a' for out- 
going transverse radiation. From a physical viewpoint, 
since (|er| 2 , |t| 2 , |t'| 2 , \a'\ 2 ) / A-kG have been identified as 
energy densities, it is perhaps more correct in general to 
encode the radiation in the spin-coefficient than the com- 
ponent of "J, since one can say that the radiation exists 
if and only if the relevant energy density is non-zero. 

It should be noted that longitudinal gravitational ra- 
diation does not exist in linearized Einstein gravity or at 
infinity, where the Bondi flux is entirely determined by 
the transverse gravitational radiation. Using the above 
expansions at future null infinity, the energy density 8n 
of outgoing transverse radiation falls off as 1/r 2 , while 
the energy density ©01 of outgoing longitudinal radia- 
tion falls off as 1/r 4 . Thus longitudinal gravitational ra- 
diation from distant astrophysical sources will be much 
too weak to measure. On the other hand, the energy 
density is generally non-zero at finite distance and does 
contribute to the irreducible mass of a black hole 0, . 

Energy conservation. The spin-components of the gen- 
erating vector £ and its dual 



(14) 



may be written as r = t^I + t 11 '!' and £ = t 00 ' I - t 11 ' I' . 
A canonical flow of time is generated by the recently in- 
troduced vector 



X = *X 



X 



-g-\dR). 



(15) 



Note potential confusion between these vectors (x, r) and 
the spin-coefficients denoted by the same symbols, but 
the distinction should be clear in context. Since D R = 
—Rp and D'R = —Rp' for uniformly expanding flows, 
one finds dR = -R(pn + p'n'), x* = R(p'l + pl')/xx an d 
X = R(p'l - pi')/ XX- Note also that x* = RH/2 and 
X = RH*/2. 

Writing the covariant derivative of the energy as 
V/:E = t 00 ' DE-t 11 ' D'E, the four bracketed terms com- 
ing from the energy propagation equations (|9"|)- (|10[l may 
all be written in a similar form: 



V f E = 



1 



4ttG 



,00', W/ft , ct 

T X l^OOO'O' + — 000'0' 



r 00 x n '($ioi'0' 



311 + Saw) 
+ r 11 'x 00 '($oio'i' +3n + S io'i') 
+ T u y i '($iii<i'+Sui'i')l • (16) 
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Collecting terms in an invariant form, 



BA'B' + 3AeAB£A'B' + ^ABA'B') 



x X AA 'r BB '. 



In tensorial notation, 



e^ a E = j>*{T ab + e ab )x a T b 



(17) 



(18) 



This has the form of an energy conservation law, express- 
ing the increase in mass-energy E due to the energy den- 
sities of matter and gravitational radiation, with respect 
to the canonical time vector \ and the normal vector t, 
scaled correspondingly to the generating vector £. It has 
the same form as the recently discovered energy conser- 
vation law for black holes Q. 

In an untrapped region where p < and p' > 0, x 
is future causal, while r is future causal if £ is outward 
achronal. Then the energy conservation law l|18|) for mat- 
ter satisfying the dominant energy condition implies the 
monotonicity property \/^E > In the special case 

of an inverse mean-curvature flow £ = ry, one has r oc x 
and the energy conservation law shows monotonicity of 
E under just the weak energy condition [l4| . 

Conclusion. An energy conservation law has been de- 
rived for uniformly expanding flows, generalizing inverse 
mean-curvature flows. It expresses the change in Hawk- 
ing energy as a surface integral of the energy densities of 
matter and gravitational radiation, with respect to a vec- 
tor x gi vm g a preferred flow of time. The gravitational 
radiation is described by an effective energy tensor O in 
the space normal to the surfaces. The Bondi energy-loss 
equation can be recovered as the limit at future null infin- 
ity, but the energy conservation law for black-hole hori- 
zons Q is independent, the expansion of (the marginal 
surfaces in) the horizon being generally non-constant on 
the marginal surfaces. Obtaining the same form of ef- 
fective energy tensor and energy conservation law is re- 
markable. 

The results hold locally for a flow of any causal type, 
generated from any initial surface on which (p,p') do 
not change sign. However, the physical interpretation 
in terms of energy cannot be so generally valid, since the 
Hawking energy in flat space-time vanishes only for met- 



ric spheres 0. The main intended application is to spa- 
tial hypersurfaces containing one or more marginal sur- 
faces. The flow gives a locally unique way to propagate a 
marginal surface outwards in the hypersurface. The level 
surfaces can presumably be found numerically by mod- 
ifying existing codes which find level surfaces of a nor- 
malized null expansion The inverse mean-curvature 
flow has been instrumental in recent proofs of the Pen- 
rose inequality in the case of a time-symmetric hypersur- 
face 0- More exactly, the flow is generalized to a weak 
flow which can jump over obstacles such as other black 
holes. If this can be implemented either analytically or 
numerically for general uniformly expanding flows in, for 
instance, a binary black-hole merger, the resulting weak 
flow could connect not only a common black-hole hori- 
zon, but also individual black-hole horizons, to infinity. 
Radiation profiles would be discontinuous at the jump, 
but such behaviour may be expected due to gravitational 
lensing by the black holes themselves. 

This energy conservation law therefore provides a mea- 
sure of gravitational radiation in the strong-field regime 
on and near a black hole, and perhaps all the way out 
to infinity. In particular, the outgoing transverse radia- 
tion has energy density 9n = |cr'| 2 /47rG. As a physical 
measure of the waveform, one may take the component 
of energy density which contributes to D' E at infinity, 
Rp<d ii /\x\ 2 , multiplied by R 2 . The only drawback of 
this proposal compared to the other made recently £| 
is that one does not have exact advanced and retarded 
coordinates, since (x,x') are defined separately at each 
moment in time. In practice, one could take appropri- 
ate linear combinations of the area radius R and the 
original time coordinate. Then the waveform-prediction 
problem reduces to whether the waveform, as a graph of 
i? 3 p6n/|x| 2 against retarded time at given angle, con- 
verges in advanced time within the numerical domain. 
Ironically, the structure required to define gravitational 
radiation is provided by the black hole itself. 
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